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Abstract 

We compute the spectrum of the su(m) spin Sutherland model of Bpj type, including the 
exact degeneracy of all energy levels. By studying the large coupling constant limit of this 
model and of its scalar counterpart, we evaluate the partition function of their associated 
spin chain of Haldane-Shastry type in closed form. With the help of the formula for the 
partition function thus obtained we study the chain's spectrum, showing that it cannot 
be obtained as a limiting case of its BCn counterpart. The structure of the partition 
function also suggests that the spectrum of the Haldane-Shastry spin chain oi Bn type 
is equivalent to that of a suitable vertex model, as is the case for its An—i counterpart, 
and that the density of its eigenvalues is normally distributed when the number of sites 
N tends to infinity. We analyze this last conjecture numerically using again the explicit 
formula for the partition function, and check its validity for several values of N and m. 

Keywords: Calogero-Sutherland spin models, Haldane-Shastry spin chains, Dunkl 
operators 

PACS: 75.10.Pq, 05.30.-d, 03.65.Fd 



1. Introduction 

The study of quantum integrable systems with dynamical degrees of freedom ex- 
hibiting long-range interactions had its origin in F. Calogero's celebrated 1971 paper (T], 
where the spectrum of an TV-particle system on the line with two-body interactions in- 
versely proportional to the square of the distance and subject to a confining harmonic 
potential was exactly computed. An exactly solvable trigonometric variant of this model 
was introduced by Sutherland soon afterwards [21 13] . The particles in this model move 
on a circle, with two-body interactions proportional to the inverse square of their chord 
distances. Both of these integrable models can be substantially generalized by exploiting 
their connection with classical root systems, uncovered by Olshanetsky and Perelomov [4]. 
More precisely, these authors noted that both the Calogero and Sutherland models are 
closely related to the -Ajv-i root system, and constructed generalizations of these models 
associated with any (extended) root system. 
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In a parallel development, Haldane and Shastry independently found an exactly solv- 
able quantum spin-| chain with long-range interactions 0(5]. The lattice sites of this 
su(2) Haldane-Shastry (HS) spin chain are equally spaced on a circle, all spins interacting 
with one another through pairwise exchange interactions inversely proportional to the 
square of their chord distances. A close relation between the HS chain with su(m) spin 
degrees of freedom and the su(ra) spin version of the Sutherland model [THS] was sub- 
sequently established using the so-called "freezing trick" (TUJ QT] . More precisely, it was 
found that in the strong coupling limit the particles in the spin Sutherland model "freeze" 
at the coordinates of the equilibrium position of the scalar part of the potential, and the 
dynamical and spin degrees of freedom decouple. The equilibrium coordinates coincide 
with the equally spaced lattice points of the HS spin chain, so that the decoupled spin 
degrees of freedom are governed by the Hamiltonian of the su(m) HS model. Moreover, 
in this freezing limit the conserved quantities of the spin Sutherland model immediately 
yield those of the HS spin chain, thereby explaining its complete integrability. The appli- 
cation of the freezing trick to the rational Calogero model with spin degrees of freedom 
also led to a new integrable spin chain with long-range interactions [10] . The sites of this 
chain — commonly known in the literature as the Polychronakos or Polychronakos-Frahm 
(PF) spin chain — are unequally spaced on a line, and in fact coincide with the zeros of 
the Hermite polynomial of degree N . The exact partition functions of both the PF 
and HS spin chains have been exactly computed by applying the freezing trick [131 114) . 

Over the years, exactly solvable and integrable one-dimensional quantum many-body 
systems with long-range interactions have attracted a great deal of attention in both 
the physics and the mathematics literature. In particular, this type of systems have 
appeared as paradigms of various condensed matter systems exhibiting generalized ex- 
clusion statistics [T5Ml7| . the quantum Hall effect [18] . and quantum electric transport 
phenomena |191 120j . More recently, quantum integrable spin chains with long-range in- 
teractions have played a key role in calculating higher-loop effects in the spectra of trace 
operators of planar JV = 4 super Yang Mills theory |2T 24 . In the mathematics litera- 
ture, this type of systems has been found to be relevant in different fields such as random 
matrix theory [35], multivariate orthogonal polynomials and Dunkl operators (2"oT429| . 
and Yangian quantum groups |30H33) . 

Spin generalizations of the BCn Calogero-Sutherland model have been extensively 
studied in the last few years, and various properties of their related spin chains of HS 
type have been analyzed with the help of the freezing trick [3~4T(4"1"] . Among the other 
classical root systems, the exceptional ones are comparatively less interesting in this 
context, since their associated models consist of at most 8 particles. On the other hand, 
although the B^, Cm and Dm scalar Calogero-Sutherland models have been studied in 
the literature [42( I43j. their spin versions have been largely ignored. Recently, however, 
the Dm spin Calogero [H] and Sutherland [25] models, as well as their associated spin 
chains, were studied by the present authors and shown to be nontrivial reductions of 
their BCn counterparts. 

More precisely, consider the Hamiltonian of the su(m) spin Sutherland model of BCn 
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type H3 Eg 

ff (BC) = _ J- d 2 Xi + a [ sin" 2 x- 3 (a - e Sy) + sin" 2 x+ (a - e5«)] 

+ b sin" 2 a; 2 (6 - e' Si) + b' ^ cos" 2 ^ (&' - e'S 1 ,) , (1) 

i i 

where the sums run from 1 to N (as always hereafter, unless otherwise stated), a, b, b' > 
1/2, e,e' = ±1, and xfj = Xi ± Xj. The operators Sij and Si in the above Hamiltonian 
act on the finite-dimensional Hilbert space 

S = ^ \si,...,s N ) | Si = -M,-M +l,...,M^, M=^-^e*, (2) 
associated with the particles' internal degrees of freedom, as follows: 

Sij \si, . . . , Si, . . . , Sj , . . . , S m) • • * j Sj ? • • • ? ^13 ' • • 3 ^n) ) 

Si\Si, . . . , Si, . . . , S N ) = \SI, . . . , -Si, . . . , S N ) , 

and we have also used the customary notation Sij = SiSjSij. If the values of the coupling 
constants b and b' in Eq. ([T]) are chosen as indicated in Table [l] one obtains su(m) spin 
Sutherland models related to the Bn, Cm and Dn root systems. 



Parameters 


Root system 


b > \ , b' = 


-Bat 


b = V>\ 


Cat 


b = b' = 





Table 1: su(m) Sutherland models of B?f, Cm, and Djv types 

Since (for instance) the Hamiltonian of the -Djv-type su(m) spin Sutherland model is 
obtained by formally setting b = b' = in its BCjq counterpart, one may naively think 
that all physically relevant properties of this -Djv-type model can also be derived from 
their corresponding BCn analogs by simply taking the (b,b r ) — > limit. However, the 
explicit computation of the spectrum of the model recently performed by the authors 
shows that this is actually not the case |45j (a similar conclusion is reached when com- 
paring the spectra of the BCm and Dm Calogero models |44j.) As a matter of fact, 
the spectrum of the -Djv-type spin Sutherland model cannot be obtained from its BCm 
counterpart |39j through any simple limiting procedure for the following two reasons. 
First of all, the Weyl-invariant extended configuration space of the Dm model — which 
turns out to be the iV-dimensional generalization of a rhombic dodecahedron — does not 
coincide with that of the BCm model, which is simply a hypercube. As a consequence, 
the (scaled) Fourier basis of the Hilbert space of the BCm model's auxiliary operator no 
longer spans a complete set of the Hilbert space of the corresponding operator for the Dm 
model. Secondly, while in the BCm case only one projector of either positive or negative 
chirality is needed to construct the Hilbert space of the model from that of its auxiliary 
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operator, two projectors of BCn type with opposite chiralities are simultaneously needed 
in order to perform a similar construction for the Dn model. Due to these two reasons, 
the Hilbert space of the Dn spin model consists of four — and not one, as in the case of 
the BCn spin model — different sectors, characterized by their chirality and parity under 
reflections of the particles' coordinates. This fundamental difference explains why the 
spectrum of the D^r-type spin Sutherland model is essentially different from that of its 
BC'n counterpart. It also accounts for the greater complexity of the partition function 
of the associated chain of Dn type [JS] compared to its BCn version [55] . 

Motivated by the nontrivial character of the -Djv-type models, one can look for other 
similar reductions of the BCV-type spin Sutherland model and their related spin chains. 
From the above remarks, it is clear that such nontrivial reductions can only be obtained 
when one or both of the parameters b and b' are equal to zero, since in that case the 
singularities of the Hamiltonian at Xi = kn and/or Xi = (2k + 1)tt/2 (with k e Z) 
disappear, so that the configuration space of the model differs from that of the general 
BCn model. In other words, the only possible nontrivial reductions of ([lj are the Dn 
model (b = b' = 0), the Bm one (b > 1/2 and b' — 0), and the model with 6 = and 
b' > 1/2. The latter model, which is not associated with a root system, is nevertheless 
equivalent to the Bn model under the change of variables Xi i— > x^ + \. Thus, apart 
from the Dn model already studied in Ref. [45] . the only new nontrivial reduction of 
the BCn Sutherland model is the Bn one. The aim of this paper is precisely to study 
the -BAr-type su(m) spin Sutherland model and its related spin chain. At the level of the 
Hamiltonians, the -BAr-type spin Sutherland model is also closely related to the Dn one, 
formally obtained by setting b = in the Bn Hamiltonian. Nevertheless, our analysis will 
reveal that the spectrum of the .BAr-type spin Sutherland model is essentially different 
from those of its -BCat and Dn counterparts. It should also be noted, in this respect, 
that the Sutherland Cat model is just a trivial reduction (i.e., a particular instance) of 
the general BCn model (fTl), obtained from it by simply setting b = b' . 

The organization of this paper is as follows. In Section [2] we define the Hamiltonians 
ff( B ) and of the spin and scalar Sutherland model of Bn type, respectively. By 

using the freezing trick, we then construct the Hamiltonian of the associated spin 
chain of HS type. We show that the sites of this chain, defined as the coordinates of the 
(unique) equilibrium point of the scalar part of the spin Hamiltonian in the principal 
Weyl alcove of the Bn root system, can be expressed in terms of the roots of a suitable 
Jacobi polynomial. Using this characterization, we establish the precise relations between 
H( B ' and the Hamiltonians of the HS spin chains associated with the BCn and Dn root 
systems. It turns out that, unlike the case of the corresponding spin Sutherland models, 
T~L^ cannothe obtained from the Hamiltonian of the -BCV-type HS spin chain by taking 
a suitable limit of its parameters. Section [3] is devoted to the computation of the spectra 
of the spin Sutherland model of Bn type and its scalar version. The main idea behind this 
computation is relating the Hamiltonians of these models to an auxiliary scalar operator 
H', which is a sum of squares of commuting Dunkl operators of Bn type. Using this 
property, we explicitly show that H' is upper triangular in the non-orthonormal basis 
introduced in Ref. |45j for the Dn model. In this way one can compute the spectrum of 
the operator H', which yields the spectra of both the scalar and spin Sutherland models 
of Bn by projecting onto suitable spaces. These results are used in Section[4]to compute 
the partition function of the HS spin chain of Bn type as the large coupling limit of 
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the quotient of the partition functions of the spin and scalar Sutherland models. Using 
the expression for the partition function of the Bjq chain, we compare its spectrum with 
those of its BCn and Dn counterparts for several values of the number of sites and 
internal degrees of freedom, verifying in this way that it is not a limiting case of the 
latter spectra. 



2. Construction of the Bjy-type HS spin chain 

Our main aim in this section is to construct the B^-type HS spin chain from its 
related spin Sutherland model by means of the freezing trick. To this end, let us first 
explicitly write down the Hamiltonian of the B^-type spin Sutherland model as 

H (B) = _J2d$.+a [ sin ~ 2 ^ (a-e S^+sin- 2 x± (o-6%)]+& ^sin" 2 ^ (b-e'S,) , 

i ijij i 

(4) 

where a,b > 1/2 and e,e' = ±1. The configuration space A^ of the Bjy model Q is 
determined by the hard-core singularities of the Hamiltonian on the hyperplanes Xi±Xj = 
kn, xi = kn (with i j and k £ Z). More precisely, we shall take A^ as the open 
subset of M. N defined by the inequalities 

< Xi ± Xj < 7r , 1 < j < i < N ; < xi < n , 1 s$ I ^ N . (5) 
It is straightforward to check that this set can be equivalently expressed as 

,4 (B) = {x e R N : < xi < x 2 < ■ ■ ■ < x N < n - x N ^} , (6) 
which is again the principal Weyl alcove of the Bm root system 

-(±e t ±e,), l^i<j^N; ±-ej, l<Z<iV. (7) 

TT TT 

Since all wavefunctions and their currents vanish on the boundaries of A^ B \ the Hamil- 
tonian if( B ) is naturally defined on some suitable dense subspace of the Hilbert space 
L 2 (A( B) ) ® S. Let us now put b = (3 a (where (3 > 0) in Eq. Q and decompose H^ B) into 
two parts as 

tf( B ) 4a fc(x), (8) 

where 

#ic B) = - E d l + «(« - 1) E [ sin ~ 2 x 7i + sin ~ 2 4 } + W a - !) E sin_2 ^ > ( 9 ) 

1 ijtj i 

which contains only dynamical degrees of freedom, is the Hamiltonian of the scalar 
Sutherland model of -Bjv-type, whereas 

M*) = \ E [ sin ~ 2 (! ~ eSii) + sin_2 4 (1 - tSii)] + f E sin_2 ^ (! - e '^) ( 10 ) 

i<j i 

is a position-dependent multiplication operator featuring the spin degrees of freedom. 
For the purposes of applying the freezing trick, we consider the a —} oo limit of ff (B ) 
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(while keeping the value of f3 fixed). The coefficient of the term of order a 2 in the r.h.s. 
of Eq. ^ , which is given by 



f/( B )(x) =^(sin" 



sin 2 x+ j ) 



sm 2 Xj 



(11) 



becomes the dominating interaction in this limit, ft was shown in Ref. |37j that the 
scalar potential t/( B )(x) has a unique minimum i9 = , . . . ,i?jv) m the configuration 
space A^ B \ By formally replacing Xi by i9i in the r.h.s. of Eq.(lO), one obtains the spin 
chain Hamiltonian 



H (B) 



sm"^r (1- eS ij )+sm~ li dJ, {1-eSi, 



J2sin-%(l-e'Si), (12) 



where df- = $i ± $j. Now, for sufficiently large a all the eigenfunctions of are 
sharply peaked around the unique minimum # of the scalar potential in the set 

@g. Hence, if ^(x) is an eigenfunction of with energy Ef c and |<jj) is an 

eigenstate of the chain Hamiltonian T~L^ with eigenvalue £j, for o>lwe have 



h(x)(pi(x)\aj) ~ v'iW^Wkj) = ^i( x )^ 



(B) 



^■^i( x )ki> 



(13) 



By using Eqs. ([8| and (13), it is easy to check that H^' is approximately diagonal in 
the basis with elements (pi(x)\crj), and its eigenvalues Eij satisfy the relation 



Ej 



ES 



4a£j 



a > 1 



(14) 



In other words, due to the decoupling of dynamical and spin degrees of freedom in the 
a — >• oo limit, the multiplication operator /i(x) can be effectively replaced by in 
Eq. pH). Consequently, in analogy with the case of other root systems, it is natural to 
define the operator T-L^ in Eq. (12) as the Hamiltonian of the HS spin chain of 2?jv 



type. At first glance, it may seem that one can use Eq. (14) to obtain each eigenvalue 



£j of the chain (12) in terms of a certain eigenvalue E^ of the spin Sutherland model 
of Bn type Q and a corresponding eigenvalue £f c of the scalar model In practice, 
however, the fact that the eigenvalues E^ and Ef c are obviously not independent makes it 



impossible to use Eq. ( 14 1 directly to determine the spectrum of the chain ( 12 ) in terms 



of the spectra of the Hamiltonians and Hsc . The key idea behind the freezing 

trick method pioneered by Polychronakos [T3] is to use Eq. (14) to compute the chain's 



partition function. Indeed, the latter equation immediately yields the exact identity 



Z{T) = lim 



Z(4aT) 
Z sc {4aT) 



(15) 



expressing the chain's partition function Z in terms of the partition functions Z and Z sc 
of the Hamiltonians and Hsc , respectively. We shall make use of this equation in 

Section [4] to explicitly compute the partition function of the HS spin chain of B^-type. 

In the rest of this section we shall discuss the relation of the B^ chain ( 12 ) with their 
BCn and Dn counterparts. To this end, recall [37] that the unique minimum of the 
scalar potential (x) in the configuration space actually coincides with the unique 
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maximum in this set of the ground state wave function of the scalar Hamiltonian 
given by 



p(x) = J^J | sincCjj sin x fj | ° | sin: 



I /3a 



(16) 



The lattice sites $j of the chain ( 12 ) arc thus the unique solution in of the nonlinear 
system of equations: 



( cot tiij + cot i?£ ) + /3 cot #i = , 1 i N . 



(17) 



In order to simplify this system, in analogy with the BCn and Djq cases let us define a 
new set of variables (fj) as 

fi = cos(2i?j) , 1 sC i sC N . 
Since # 6 A^ B \ from Eq. ^ we obviously have 



< i?i < •• • < i?iv_i < 



< $ N < 7T 



and therefore 



1 > fi > 6 > • • • > Civ-i > , f w < i 



In terms of the variables (fj), the system (17) can be rewritten as 



6-6 1-6 



= 0. 



1 sC i < N. 



(18) 



(19) 



Let fj = min{fi, . . . ,£ N } = min{fjv-i, £ N }. Since 1 - f io > and fj - f io > for all 
j 7^ «0) the system (19 1 with i = io implies that fj = — 1. From Eq. (18) it follows that 
iq = N. Substituting fjy = —1 into (19), we obtain the following system of equations for 
the remaining coordinates fi, fa, . . . , £n— i' 

N-l 

2 ( 1 -^E7T7 = ^ 2 + ^ + 2 )' l<t<JV-l. (20) 



Comparing ( 20 ) with the system 

JV 



2(1 - C?) E T^TT =P-P + (P + P% , 1 < » < , 



(21) 



satisfied by the zeros Ci (i = 1, . . ■ , N) of the Jacobi polynomial hP ' 1] (cf. Ref. @7]), 

we conclude that the coordinates fi,f2, ■ ■ ■ ,fiv-i are the zeros of P^Zx ■ In terms of 
the original site coordinates i?j we have 



0<#1 <tf 2 < ••• <0N-1 



N 



with F \ 



08-1,1), 



N-l 



3(2i?i)) = for % = 1,2,... ,N- 1. 



For the purpose of comparing the B/v-type HS Hamiltonian ( 12 ) with its BCn coun 



terpart, let us now briefly review the construction of the latter spin chain from the 
corresponding spin dynamical model 39J. Due to the singularities at the hyperplanes 



Xi±Xj = kn, Xi = kn and Xi 



- kn (with 1 ^ i < j ^ N and k G Z) , the configuration 



space of the spin Sutherland Hamiltonian ([I]) can be taken as the principal Weyl alcove 
of the BCn root system 



^(BC) 



< Xi < X2 < 



< xn < 



'■}■ 



(22) 



Applying the freezing trick to the Hamiltonian ([!]) with b — f3a and b' — j3'a (with 
> 0), one obtains the su(m) HS spin chain of BCn type as 



^ BC) = *£[sirr 2 0r.(l- e ^ 

i<j 



-2 a+ 



4 2-j 



cos 



■Oiil-e'Si), (23) 



where 9^ = w % _l_ v 3 
the scalar potential 



± 0j and = (6*1, ... , is the unique equilibrium in the set A^ BC ^ of 



C/( BC >(x)=£(si 



2 - 

sm x„- 



sin 2 x+-) + >(/3 2 sin~ 



: + P' 2 COS 2 Xi) 



(24) 



In fact, it is shown in Ref. [37] that the chain sites 9^ can be expressed as Q = cos(2#j), 

where Q are the zeros of the Jacobi polynomial P^ 1,13 . 

Let us now try to find out the precise relation between the .B/v-type HS spin chain 



(12) and the j3' — > limit of its BCn counterpart (23). In this context, it should be 



noted that the potential £/( BC ) in Eq. (24) smoothly reduces to the Bn potential 



in Eq. (11) in the limit /?' — > 0. Consequently, the lattice points of the spin chain (23) 



should coincide with those of the spin chain (12) in the /3' — > limit. In other words 



Eq. ( 19 1 should also yield an alternative characterization of the coordinates as the 



TV roots of the Jacobi polynomial P$ l ' . This is indeed the case, since Eq. (19) 



obviously reduces to (21) when j3' — 0. Alternatively, by using well-known properties of 



the Jacobi polynomials we can easily establish the identity 



P 



N 



(3 + N- 



2N 



(25) 



which confirms the equivalence of both characterizations of the site coordinates of the 
■Bjv-type spin chain. 

Next, using the identity cos -2 9{ = 2/(1 + d), we can express the last term in the 
r.h.s of Eq. pi as 



1 



2 4- i + d 



(1 - e'Si) . 



(26) 



Since Q — > £j as j3' —> 0, it is clear that the relation lim^'^oCl + d) — holds only for 
i = N. As a result, all the terms but the last one in the sum in Eq. (26) tend to zero as 



j3' — > 0. In order to evaluate the limit of this last term, we divide (21 I by 1 + Q and sum 
the resulting equation over i, obtaining 



•c 



Taking the f3' — > limit of both sides of the above equation, and using the fact that 
Bmp^o P'/(l + Q) = for i ^ AT, we get 



lim 



2/3' 

roiTc" 



N 



N(N-l + 0). 



From Eqs. (23), (26) and (27) it immediately follows that 



lim n (BC) = H (B) + - N(N - 1 + p) (1 - 6'Sat) . 



(27) 



(28) 



Thus the /?' — > of the Hamiltonian of the HS chain of BC^ type yields its analog, 
plus an additional term which can be interpreted as an "impurity" at the right end of 
the latter chain. 

Consider now the Hamiltonian of the su(m) spin Sutherland model of Dn type, which 
is obtained by setting b = b' = in Eq. ([lj: 

ff (D) = _ J-dl + a [ sin- 2 z7. (a - e S tj ) + sin" 2 x% (a - e S tf )] . (29) 

The configuration space A^ D ^ of the model (29) is determined by the hard-core 
singularities of the Hamiltonian on the hyperplanes Xi ±Xj = kir (with i ^ j and k € Z). 



For TV > 2, it is easy to check that is given by 

A (D) = {x e R N : \xi\ < x 2 < ■ ■ ■ < x N < 7r - x N -i} . 



(30) 



which is again the principal Weyl alcove of the Dm root system. Application of the 
freezing trick to the Hamiltonian ( 29 ) leads to the Hamiltonian of the su(m) HS spin 
chain of type given by 



H <P) 



i<j 



(l-eSy)+sin- 2 0^(l-aS, 



(31) 



where the lattice sites & are the coordinates of the unique minimum cj) in the set 
of the scalar potential 



C/( D )(x) = ^(sin- 

»*7 



sin 2 x+ 



(32) 



Again, defining new variables Xi — cos(2^), one can show that xi — ~Xn = 1 and that 
the coordinates X2> ■•■ iXN-i are the zeros of the Jacobi polynomial P^l- Using this 
characterization, it has been shown that in the (/?,/?') — > limit the Hamiltonian of the 
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HS chain of BCn type yields its Dn 
latter chain 1451: 



analog, plus "impurity" terms at both ends of the 



hm H^=H^ + -N(N-1) 
WHO 2 



-(Si + Sn) 



(33) 



Let us now try to establish a relation between the Hamiltonians of the HS spin chains of 



-B/v and Dn type. To this end, we take the /? — > limit of both sides of Eq. (28), which 
yields 

lim H {BG) = lim U {B) + ] N(N - 1) (1 - e'S N ) . 



WHO 



/3^0 



Comparing the r.h.s. of the above equation with that of Eq. (33 1 we obtain the relation 

(34) 



lim H {B) = H (D) + \ N(N - 1) (1 - e'Si) , 

p^o 4 



which shows that the /3 — > limit of the Hamiltonian of the HS chain of Bn type yields 
its Dn analog, plus an "impurity" term at the left end of the latter chain. 

As mentioned earlier, the Hamiltonians of the Bn and Dn spin Sutherland models 
can be obtained from their BCn counterpart by formally taking some limits of the 
related coupling constants. On the other hand, due to the presence of impurity terms 
in Eqs. (28 1, (33) and (34), it is clear that the Hamiltonians of the HS spin chains 
associated with the Bn, Dn and BCn root systems cannot be related to each other by 
any simple limiting procedure. Hence, it is natural to expect that the spectrum of the 
HS spin chain of Bn type should be qualitatively different from both its BCn and Dn 
counterparts. In this context it should be noted that, in spite of the apparent closeness at 
the level of their Hamiltonians, the configuration spaces of the Bn, Dn and -BCV-type 
spin Sutherland models are completely different from each other. Indeed, by comparing 
Eqs. ([§, ([22]) and ((30| with each other, we find that D D A^ BC \ Since the 



Hilbert space of a dynamical model is built up from square-integrable functions defined 
on the corresponding configuration space, this result clearly indicates that the spectrum 
of the spin Sutherland model of Bn type should be qualitatively different from those of 
both its BCn and Dn counterparts. 



3. Spectra of the Bjv-type spin Sutherland model and its scalar version 

In this section we shall compute the spectra of the spin Sutherland model of Bn 
type |3and its associated scalar model We shall employ a well-established tech- 
nique [301 [Ml 55], which is based on relating both of these models to the auxiliary 
differential-difference operator 

h' = - ^2 d li + a X] [ sin ~ 2 x ij (° ~ K v ) + sin ~ 2 x tj ( a ~ ^« )] + b s[n ~ 2x i 

i t/j i 

(35) 

where and Ki are coordinate permutation and sign reversing operators, defined by 



(■Kijf) (^Xi ■ • ■ ■ ■ • i^ji • • • , ^N) — fy^Xi • ■ • , ■ • ■ ■ 

(K l f)(xi, ...,Xi,..., x N ) = f(xi, -Xi, x N ) , 
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,xn) , 



and Kij = KiKjKij. Let us denote by 227 the group generated by the operators Km 
and Ki, i.e., the Weyl group of the BCn root system, which actually coincides with that 



of the -Bjv and Cm systems. From Eq. (35) it is clear that the operator H' is naturally 
defined on a dense subset of L 2 (C (B) ), where C (B) = 2B • A^ B K We shall next show that 

C (B) = {x G R N : < \xi ± Xj \ < tt, x 4 ^0; 1 s? i < j s? N} . (36) 

Indeed, first of all it is obvious that 

217- A (B) C {x G R N : < |x; ±xA < tt , 0<|x;|<7r; 1 ^ i < j ^ N} . 

Adding the two inequalities — tt < x, ± Xj < tt we immediately obtain — ir < Xi < w, so 
that 

2U ■ A {B) c {x G R N : < |x* ± xj\ < tt , x t ^ ; 1 s? i < j s? TV} = C* . 

Hence to prove p6| we need only show that C C W-A^ B \ To this end, note that if x G C 
there is an element W of the Weyl group 2U such that Wx = y, where < y\ < ■ ■ ■ < y^. 
Since C is invariant under 2U, the vector y must belong to C, so that yw-i + Vn < 7r - 
Hence y G A (B) , and therefore x = W^y G 2U • A^ B \ This shows that C C 2U • A( B \ 



thus completing the proof of Eq. ( 36 ) 



As mentioned in the previous section, the operators and are naturally 

defined on suitable dense subspaces of the Hilbert spaces L 2 (A^ B ^) <g> £ and L 2 (A^), 
respectively. In order to compute the spectra of these operators, we shall start by con- 
structing suitable isospectral extensions H and H sc thereof to appropriate subspaces of 
L 2 (C (B) ) ® £ and L 2 (C (B) ), such that H = H 1 ®1 and # sc = if' in the latter subspaces. 
More precisely, denote by A ee > the projector onto the subspace of L 2 {C^)®H consisting 
of states with parities e and e' under particle permutations and simultaneous reversal 
of each particle's coordinate and spin, respectively. In other words, the operator A ee < 
satisfies the relations 

KijSijA^, = ek ee , , KiSik eel = e'K ee , . (37) 

As shown in Ref. [35], there is a natural isomorphism ~ between the spaces L 2 (A^ B ^)(g)S 
and A ee / (L 2 (C (B) )(g)S) , so that iJ (B) is isospectral with the operator H = ~oij( B )o(~)- 1 
defined on an appropriate dense subset of the latter space. Similarly, if we denote by A sc 
the projector from L 2 {C^) onto the space of functions symmetric under permutations 
and even under sign reversals, defined by the relations 

KijKc = KiA sc = A sc , (38) 

the spaces L 2 (^4( B )) and A sc (i 2 (C ( - B ^)) are again naturally isomorphic. Hence, denoting 
(with a slight abuse of notation) this isomorphism by ~ , the operators H sc and H SI 



o H sc o ( ) 1 are again isospectral. From Eqs. ( 37 )-( 38 ) and the definition (35) of the 
auxiliary operator H' , it immediately follows that 

H = H' ® l| Aee ,( i 2( C .(B) )!glS ) i Hsc = #'Ia sc (l2( C (b))) ■ ( 39 ) 

In order to compute the spectra of H and H sc , we shall first triangularize the auxiliary 
operator H' , whose domain is (a dense subset of) the Hilbert space L 2 {C^). In fact, 
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L 2 (C^) ee L 2 (C( B )), where denotes the closure of the set C {B \ Using Eq. @ 

it is immediate to show that 

C (B) = {x G R N : \xi ± ar^l < tt, 1 s$ i < j < TV} , 

which coincides with the analogous set for the spin Sutherland model of Dn type studied 
in Ref. [45 . As shown in the latter reference, this set is the TV-dimensional version of a 
rhombic dodecahedron. Furthermore, it was shown in Ref. (45) that one can construct 
a basis of the Hilbert space L 2 (C^ B ^) by considering the complex exponentials e 1 ^ kjXj 
(with (ki, .. . , k^) G Z N ) which are periodic on C^ B \ namely the set of functions 

e iE,(a«j+i)»i ) n = (m,...,n N ) eZ, S G {0, 1} . (40) 

By using standard arguments, it can be readily proved that the set of "gauged" Fourier 
functions 

^(x)EE P (x)e i ^^+^, n=(n 1 ,...,n JV )eZ, S G {0, 1} , (41) 
where p is defined in Eq. (16), is a (non-orthogonal) basis of L 2 (C^ S '). 

5. i. Triangularization of H' 

We shall next define a suitable order in the set (41) so that the action of H' on the 
resulting basis is triangular. Note, first of all, that 

L 2 (CW)=SjW®%W, (42) 

where S)^ is the closure of the subspace spanned by the basis functions ifin^ with n G 7L N . 
We will show that H' leaves invariant each of the subspaces & s \ so that we need only 
order each subbasis {tfn } neZ w m such a way that H' is represented by a triangular 
matrix in $ry- s >. To this end, given a multiindex p ee (p 1; . . . ,pn) G Z n wc define 

[P] = (\Ph !>•••> \Pi N \) > with IPiil > ■■■ > \Pw\- 

If p' G 1> N is another multiindex, we shall write p -< p' provided that the first non- 
vanishing component of [p'] — [p] is positive. The basis functions f } nGZ jv should 

then be ordered in any way such that precedes ip^) whenever v -< u 1 , where 



v ee (2m + 6, . . . , 2n N + S) , (43) 

and similarly for v' . For instance, ) must precede <P^2 > _ 3 _^ and \ ) should fol- 
low ( p^_ 3 -1)' wnue the relative precedence of <P(2^-3 -1) an< ^ V 13 -2) can ^ e arbitrarily 
assigned. 

In order to compute the action of H' on the basis functions (41 ), we shall express the 
latter operator in terms of the Dunkl operators of B n type 



Jk = id Xk + [(! - [cotx ki) K ki + (! - ' lcotx ti) K kl 

l^k 

-2a^K kl +b{l-icotx k )K kl k = l,...,N, (44) 



Kk 
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obtained from their BCn counterparts in Ref. [39] by setting b' = 0. Note that the 
natural domain of the operators J k is the same as that of H', i.e., a suitable dense 

=?(B) 



subspace of L 2 (& '). Setting b' = in Eq. (10) of Ref. [39] we obtain 



(45) 



so that the action of H' on the basis ( 41 ) can be deduced from that of the Dunkl 
operators (44). In the following discussion, we shall label the basis functions ipn simply 
by tp u , with v defined by (43). As in Ref. [15] . we shall start by considering the action of 
Jk on a basis function <p v with v nonnegative and nonincreasing. For such a multiindcx, 
we shall use the notation 

#(s) = card{i : Vi = s} , £(s) — min{i : Vi = s] , 

with £(s) = +oo if Vi ^ s for all i = 1, . . . , N. For instance, if v = (9, 7, 5, 5, 3, 3) then 
#(5) = 2 and ^(5) = 3. 

Our next step is to prove the key formula 



E 



c u,k Vw 



(46) 



v' -ve{2Z) N , v'< v 



where c£ k £ C and 



^v,k — 



-v k + 2a(2t{v k ) + #K) - k - N - 1) - b, v k >0 
2a{N -k)+b, v k =Q. 



(47) 



Indeed, using Eq. (44), and performing a lengthy but otherwise straightforward calcula- 
tion one finds that 



= -v k - 2a(N - 1) + 2a^2 



a /k - 1 



j<k a i k 



2«E 

j>k 



*jk 



nl—Uj—U h 

2aJ2 - +2b Z -± 



2{l-u k ) 



3jk 



-b, (48) 



where 



otik = z i z k 



Pj k z j z k 7 



Note that all the terms in Eq. ( 48 1 except for the last two also appear in the corresponding 
formula for the Dff case, cf. |45[ Eq. (51)]. Since 



<Pi 



m-v h ) _ , 

z k L 



-<PvZ k 



2(l-w fc ) 1 ~ z k 



(z fc 2 + ■ ■ • + z. 



-2(t/ fe -l) 



Vk = 

)(p v -<<p v , u k j^0, 



(49) 
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the contribution to X v k of the terms proportional to b in Eq. (48) is equal to b(25 Vk q — 1). 
Taking this into account, together with Eqs. (49)-(50) of Ref. |45j for the case, we 
easily obtain Eqs. p6|-|27|. 

Since Eq. ( 46 ) does not hold in general when u does not belong to [Z N ~\ , Eq. ( 47 ) 
does not give the complete spectrum of the Dunkl operators J k . However, in order to 
compute the spectrum of H' we shall only need the following weaker result: if v e 7L N is 
a multiindex all of whose components have the same parity, then 



E 



(50) 



for some complex constants 7^ k . In order to prove this formula, note that if v is as 
stated above there is an element W € 2D such that <p v — W<p\ v i ■ Setting b' — in the 
commutation relations between the -BCjy-type Dunkl operators and the generators of 2D 
listed in Ref. [SS] , it is straightforward to show that 

2 N JV! 

[J k ,W] = c ikWj, c jk eR, 

where we have denoted by Wj (with j = 1, . . . , 2 N N\) an arbitrary element of 2D. From 
the previous equation and the relation ip v = WtpuA we easily obtain 

2 N N\ 

J k tp v = W{Jkf[ u ]) + ^ CjkWjPlv] ■ 
i=i 



Applying Eq. ( 46 ) to the multiindex [v] , and taking into account that the partial ordering 
-< and the parity of the components are invariant under the action of 2D, we easily arrive 
at Eq. ((501). 



We shall next show that the action of H' on each subbasis {v?n } ngzJV , ordered as 
explained above, is upper triangular: 



H'^ = BWpff) + £ <#> n¥ #> , v h ee 2n k + S , v' k ee 2n' k + 5 , 



(51) 



V 1 



where cf}_ £ C and 



(52) 



Indeed, suppose first that the multiindex v in Eq. (51 ) is nonnegative and nonincreasing 
App 
that 



Applying J k to both sides of Eq. (|46| and using Eq. ([50]) , it is straightforward to show 



J k Vv = XtfiVv + E h »,k Pi 



U < V 
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with 6£ k £ C. From the identity ( [45] ) we thus obtain 



E 

-vG(2Z) 



(E<0 

fe 



(53) 



Suppose, next, that ^ and let again € 2B be such that tp v = Wip[„y As shown 

in Ref. [39], the -BCV counterpart of the operator H' commutes with all the elements 
of 2U. Since H' is obtained from its BCm analog by setting b' = 0, it follows that 
[H', W] = 0. Using this fact and applying Eq. |53| ) to tp^ we find that 



H'<p v = W-H , <p [lf] = (E A M,*V 



E 

i/'-Me(2Z)' 



which establishes (51) with 



E n = E A fUfe 



(54) 



The last step in the proof of Eqs. (51)-(52) is to show that Eq. (54) can be simplified 



to yield Eq. (52). For this purpose, let us write p = [v] and consider first the case in 
which p k ^i > pk = ■■■ = Pk+r > Vk+r+i > 0. Since £(pk+j) = k and #(p k+J ) = r + 1 for 
0,. 



j = U, . . . , r, using Eq. (47) we obtain 
A p ,/c+j = —pk+j + 2a(k+r-j-N)-b = -pk+r—j + 2a(k + r — j - N) 



3 = 0, • 



Thus in this case wc have 



k -\-r fe+r 

E A P,; = Efe+ 2fl ( Ar --?) + & ) i 

l—k j—k 



(55) 



On the other hand, for the case in which pt-i > Pk = • • • = Pn = the analog of Eq. ( 55 ) 
follows directly from Eq. (47). Thus Eq. (55) is valid for arbitrary u £ 1 N , and Eq. (52) 



follows from the latter equation by summing over k. 

3.2. Triangularization of and 

We shall next make use of the previous results to triangularize ij( B ) and i4 B) . As 
mentioned above, this problem is equivalent to the triangularization of the extensions H 

and H sc acting on their respective Hilbert spaces S) = A ee > [L 2 (C ) <8> SJ and Sj sc = 



(B) 

A sc (i 2 (C ))j which can be carried out without difficulty with the help of Eq. (39). 

Let us start with the operator H. By Eq. ( [42] ) , its Hilbert space can be decomposed 
as the direct sum 

Si = A £e , (SjW ® E) A ££ , (i5 (1) ® E) . (56) 

Let /(x) be a function in the domain of H', and let \s) £ E denote an arbitrary spin 
state. Since H coincides with H' <S> 1 on Si, and the latter operator commutes with A ee < 
(indeed, it commutes with all the elements of 22J) , we have 



#[A ££ ,(/(x)| s ))] =A ££ ,[(#'/(x))| S >]. 
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(57) 



As H' preserves each subspace the latter equation implies that both subspaces 

^-££'(^3 ® S) with 6 = 0, 1 are invariant under i?. We shall next verify that H acts 
triangularly on a (non-orthogonal) basis of A ee < (fj'-' 5 * 1 <8 S) of the form 

<)(x)=A ££ ,(^)(x)|s)), (58) 
ordered in such a way that -0„*s precedes Vv s' whenever v < v 1 (with v defined in (431, 



and similarly u'). Since the functions tp^ are a basis of L 2 (& '), the spin wavefunc 



(B) 



tions (58) are obviously a complete set, but they will not be linearly independent unless 



suitable restrictions on the quantum numbers (n, s) are imposed. More precisely, the 
states (58) are a (non-orthogonal) basis of the Hilbert space A EE i($j^ ® S) provided 
that n £ Z w and s £ {— Af, —A/ + 1. . . . , Af}^ satisfy the following conditions: 



i) n.i ^ • • • ^ riAr ^ 

ii) If 6 = rii = then ^ for e' = 1, while Si > for e' 



iii) If Uj = rij and i < j then 



$i ^ S-j 



Si > s 



for e = 1 
for e = — 1 



(59a) 
(59b) 

(59c) 



(In condition (ii), it is understood that no additional restriction is imposed on Sj when 
either 5 or rii is nonzero). 
Indeed, since 



A ee ,{K l3 S tj ) = eA e 



A eE ,{K i S i ) = s'A ee , , 



acting with suitable operators KiSi and KijSij on a spin function tp„ i (x)|s) with arbi- 
trary n £ Z w and s one can easily show that the corresponding state ipn}s is either zero or 



proportional to a state (58) satisfying the above conditions. (Note, in this respect, that 
a state (58) with 6 = rii = Si = is symmetric under (xi,Si) — > {— Xi, — Sj), and must 
therefore vanish identically if e' = — 1.) This shows that the states (58) with n £ Z N 
and s satisfying the conditions (59) are complete. Their linear independence is easily 
checked. 

Let us now examine the action of the operator H on the basis of A SE > (fy^ ® S) given 



by Eqs. (|58|-(|59|). From Eqs. (|51])-(|52j) and the identity ( |57[ ) one immediately obtains 

(60) 



where c„) _ £ 



and 



(61) 



Although the quantum numbers (n',s) appearing in the r.h.s of Eq. (60) need not satisfy 



conditions ( |59| ), there is an element W £ 2B such that (Wn',Ws) = (n",s") do satisfy 

ling 
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the latter conditions. Since the corresponding state ip^,} s „ differs from ip^, B by at most 



an overall sign, and [v"] = [u'\ -< [v] implies that v" -< u, it is clear that we can 



rewrite ( 60 ) in the form 



rn,s n,srn,s 



y c^), 

/ j n s ,ns 



(62) 



n' s' ,ns 



Hence the action of H on each subbasis ( 58 )-( 59 ) 



for suitable complex constants c, 
(with fixed d G {0,1}) is indeed triangular, with eigenvalues given by Eq. (61 1. The 



spectrum of H is thus obtained from Eq. (61 ) when 6 = 0,1 and (n, s) are any quantum 



numbers satisfying conditions ( 59 ) . 

Since, as mentioned at the beginning of this section, the operator is isospectral 



to H, Eq. (61 1 gives the complete spectrum of the spin Sutherland model of Bn type. In 



particular, the energies of this model do not depend on the quantum number s. Therefore, 



the degeneracy dn of the eigenvalue (|61|) due to the spin degrees of freedom is simply the 



number of basic spin states |s) satisfying conditions (59). In order to explicitly compute 



this degeneracy, let us write the quantum number n in the form 

n= (pi,...,pi,...,p r ,...,p r ), 

Using conditions (59b])- ([59c]) we easily find that 

r-i 



Pi > ■ ■ ■ > p r ^ . 



d^ 



m ee >(k r ) 



n 



m £ (ki 



n 



m £ (ki 
k; 



Pr 



otherwise, 



(63) 



(64) 



where 



m e (ki) = m+ -(1 + e)(k l 



m EE ,(k r ) = - [m + e'n(m) + (1 + e)(k r - 1)] (65) 



and 7r(m) = m (mod 2) is the parity of m. 

Similarly, the spectrum of the scalar Hamiltonian H sc can be computed using the fact 

that it coincides with H' in the Hilbert space $) sc = A sc (L 2 (C ) ) , which by Eq. ( 42 ) 
is given by 

-Ssc = A sc (^ >)©A sc (i3«). (66) 

Due to the identity 

# sc (A sc /(x)) = A sc (ff7(x)) , 

it is immediate to show that each of the subspaces A sc (f)(' 5 " ) ) is invariant under ^ sc . Just 
as for the spin model (cf. Eq. (58)), it can be verified that the functions 



^(x)=A sc (^ ] (x)) 



where n € 7L N and 



ni ^ 



> n N > , 
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(67) 
(68) 



are a (non-orthogonal) basis of A sc feW). Proceeding as above, it is straightforward 
to show that the operator H sc acts triangularly on the basis obtained by ordering the 
set ( 67 1-( 68 ) in such a way that ipn^ precedes whenever v < u' , and that its 

eigenvalues E^ are given by the r.h.s of Eq. (61 ). Since the operators and H sc are 
isospectral, Eq. (61 ) gives the full spectrum of the scalar Sutherland model of Bjq type. 
Of course, due to the absence of internal degrees of freedom, in this case the degeneracy 
factors dn are equal to one for all quantum numbers n and 8 = 0, 1. 

A remark about the differences between the Hilbert space and spectrum of the spin 
Sutherland model of type and those of its BCn and counterparts is in order at 
this point. In the first place, the Hilbert space of the BC^ spin Sutherland model is 
isomorphic to the subspace A ee < ($j(°> <S) S) of Sj. On the other hand, the parameter e' is 
not present in the Hamiltonian of the Dn model, but is instead a quantum number which 
can take the two values ±1. For this reason, the Hilbert space of the spin Sutherland 
model of Dm type consists of four sectors, each of which is isomorphic to a space of the 
form K ee ,{^ ® S) with £ = 0,1 and e' = ±1. Stated differently, in the BC N model 
both parameters 8 (= 0) and e' are fixed (and, therefore, do not appear in the spectrum 
as quantum numbers), in the Bn model 8 is a quantum number but e' is fixed by the 
Hamiltonian, whereas in the Djq model both 8 and e' appear as quantum numbers in 
the spectrum. 



4. Partition function of the B^-type HS spin chain 



The purpose of this section is to evaluate in closed form the partition function of the 
Haldane-Shastry spin chain of Bn type (12) using the freezing trick. To this end, we 



shall make use of the key relation (151 expressing the chain's partition function Z in 



terms of the partition functions Z and Z sc of the Hamiltonians H ( B ) and -fficf ' . In order 
to compute the a — > oo limits of Z(4aT) and Z sc (4aT), we start by expanding Eq. (61) 



for the energies of both and Hsc in powers of a, with the result 

Egl = a 2 E + 8aJ2 n k(N + ^ - fc) + 2a SN{N + f3 - 1) + 0(1) . 



(69) 



where 



E = N/3 2 + 2N(N -l)f3+ - N(N - 1)(2N - 1) . 



Note that, since a 2 Eg does not depend on rik, it will clearly not contribute to the quotient 
Z (AaT) / Z sc (4:aT) . We can therefore subtract this term from the spectra of both H^ B > 



and -ff sc for the purposes of computing Z through Eq. (15). With this normalization 



the eigenvalues of and become O(a) for a —> oo, so that the limits of Z(4aT) 
and Z sc (AaT) exist separately. Dropping the term a 2 E from Eq. (69) we thus obtain 



lim Z{4aT) 



E 



E 



i<SJV(iV-W?-l)+2 £ miN+Z-i) 

c4 4 > q 



q = e 



5=0,1 ni^-^njv^O 



(70) 
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Using Eq. ( 63 ) it can be easily shown that 



N (3 r 

^2m{N+ - - ij = v x 



1=1 



where we have defined 



(71) 



(72) 



Substituting Eq. (|71j) in Eq. (JTOJ) we have 

limZ(4aT)=y V dg» q£™ + qi W+P-V V V dg> 

= Z (q) + Zx(q) , (73) 



where "Pat is the set of all partitions of the integer N taking order into account, and we 
have denoted by Zo(q) (resp. Z\{qj) the contribution of the 5 = (resp. 6 = 1) sector to 
lim^oo Z{AaT). 

We shall next proceed to simplify each of the functions Z& (q) . In the first place, using 
the definition of Z (q) and the value of the degeneracy factors d^ in Eq. (64) we easily 
arrive at the formula 



_ , . , , \ - 2 T, Pi "J r • \ - 2 S Pin 

z (q) = }^ d d 1 2^ q 1=1 +2^ do<fc 2^ « !=1 

kev N pi>->p,>o kePjv pi>--->p r _i>o 



(74) 



where we have set 

r-1 



i=l 



m e (ki 
k; 



d x = 



m e (k r ) 



d 2 = 



Proceeding as in Ref. [39] it is straightforward to obtain the key identity 

s 

n i q£(Ki< ■ 



E ? ,; 

Pi>— >p»>0 



2 £ Pfi A 



i=l 



where the dispersion relation is defined by 

£{t) =t(2N + 13- 1-t). 



Substituting this identity with s = r and s = r — 1 in Eq. ( 74 ) we find that 
kePjv «=i ke-Pj\ 



z o(q) = E d ° dl n i 



e n i 



(75) 



(76) 



(77) 
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Consider now the function Z\(q), explicitly given by 

This formula can be simplified by using the identity 

2Ek»i -r-r (?' 



E ^ «= 

Pl> — >Pr>0 



S{Ki 



i=l 



1 



1 - q £(Kr) 



which easily follows from Eq. ( 75 ) , with the result 

J — 1 £(Ki 

k - r 

Note that K r = N, so that 



Z 1 (q) = q ^ N +^ £ dodillr 



i=l 



_ q£(Ki) 1 _ g£(Kr) 



(78) 



£(X r ) = iV(iV + /3 — 1) 



Substituting Eqs. (77J and (78j in Eq. (|73j) we obtain 

kePiv i=l kGPjv i=l ^ 



±JV(iV+/3- 



r-1 



x) e n i 



' 1 q £(Ki) 



1 



7 £(ifi) 1 _ „£(K r ) 



e 



kePjy i=l i=l 



r 1 

II 



• do 



d 1 (gl^)+^))+d 2 (l-/(^)) 



After a straightforward simplification, this equation yields the following explicit formula 
for the a — > oo limit of the partition function of the spin Sutherland model of B^ type: 

lim Z(4aT) 

a— >oo 

= (l+ g fW-D) £ H Q £iKi) -fl x _ \ {Ki) ■ do \d iq ^ N +e-V+d 2 (l- q f^-») 



kev N i=i 



»=i 



(79) 



We shall next evaluate the partition function of the scalar Sutherland model of Bn 
type in the limit a — > oo. As mentioned above, the energies of this model are still given 
by the r.h.s. of Eq. (69), although in this case there is no degeneracy due to the spin 
degrees of freedom. Thus the large a limit of the partition function Z sc (4aT) is given by 

Eq. ([70} with 4f } = 1: 

lim Z sc (4aT) = V V q ^^g^ N+ ^ 

(I — ^-oo ' ^ ^ 

(5=0,1 rii^-^njv^O 



+ E « 



2 £ n«(JV+#-t) 



(80) 
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Evaluating the last sum as in Ref. [39] we readily obtain 



N 



lim Z sc (4aT) = (l + q 



i=l 



(81) 



The partition function of the Haldane-Shastry spin chain of Bn type ( 12 ) is easily 



computed by inserting Eqs. (79) and (81 1 into the key relation (15). In order to sim- 
plify the resulting expression, we define N — r integers K[ < ■ ■ ■ < K' N _ r in the range 
1 V 1 by 



{K{,...,K' N _ r } = {l,...,N-l}-{K 1 ,...,K r -i}. 
Using this notation, we finally arrive at the following closed-form expression for the 



partition function of the spin chain (12): 



r-1 



N — 7 



E 

kePjv i=l 



l[ q £m.l[(i- q m)).do 



diq 



f (JV+^-l) 



do 1 



(82) 



In particular, from the latter equation it is clear that Z is a finite sum of powers of q, as 
should be the case for a finite system. 

For comparison purposes, we note that the partition function Z^ BC ^ of the HS spin 
chain of BCn type (23), which can be inferred from Eqs. (52)-(54) in Ref. [US], may be 
written in the notation of the present paper as 



£(BC) 



E 



i=l 



N-r 

n 



1-* 



d x q 



N(N+P+P'-1) _ 



do 1- 



N(N+P+P'-1)\ 



(83) 

with £ (t) = t(2N + ft + /3' -1-t). Comparing Eqs. ([82]) and ([83]), it is apparent that the 
partition function ij( BC ) does not tend to its Bjq analog Z in the limit /3' — > 0. Likewise, 
it is clear that the expression of the partition function of the spin Sutherland model of 
Dn type given by Eqs. (92) and (95) in Ref. [45] is much more complex in nature than 
its Bn counterpart (82) with j3 — 0. Indeed, the fact that the partition functions of the 
BCn, Bn and models cannot be obtained from one another by taking appropriate 
limits of the parameters (3 and /3' is in agreement with the presence of boundary terms in 
Eqs. (33)-(34). In order to illustrate this remark, in Fig.[T]we have compared the spectra 
of the Bn chain with its BCn and Dn counterparts for different choices of N, m, and 
p. More precisely, in the latter figure we have plotted the (normalized) cumulative level 
density of these chains, defined by 



F(E) = 



1 



Ei^E 



where E\ < ■ ■ ■ < E n are the distinct energies and Si denotes the degeneracy of the 
energy Ei. It is apparent from these and similar plots that the spectra of the Bn, BCn 
and Dn chains cannot be obtained from one another by taking appropriate limits of the 
parameters j3 and /?'. 
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F(E) F(E) 




E 



Figure 1: Left: cumulative level density of the ferromagnetic HS chain of Bjv type with N = 12, m = 2, 
and f} = 2 (blue) vs. its BCpf counterpart with + /3' = 2 (red). Right: cumulative level density of the 
ferromagnetic HS chain of -Djv type with iV = 10 and m = 3 (red) compared to its _Bjy analogs with 
P — > and e' = 1 (blue), e' = — 1 (green). Note that, by Eq. J65J, the energies of both the Bjv and the 
BCV chains are independent of e' when m is even, while the spectrum of the -BCjv chain depends on (3 
and /3' through the combination /3 + /3' on account of Eq. (|83| . 



On the other hand, the obvious structural similarity between Eqs. (82 1 and (83) and 
the faclQ that the spectrum of the HS spin chain of BCn type can be described in terms 
of a suitable generalization of Haldane's motifs [?5] suggests that a similar description 
should also exist for the present chain. Note that, for HS chains of Am type, the existence 
of such a description is the key ingredient in the proof of the Gaussian character of their 
level density |49j when the number of sites tends to infinity, which is of importance in 
the context of quantum chaos and in the study of the thermodynamic properties of these 
chains [50 . In fact, using Eq. (82) we have numerically checked that the level density 
of the HS chain of type is approximately Gaussian when N > 10, for a wide range 
of values of the parameter /? and the spin degrees of freedom m (see, e.g., Fig. [2]). This 
property of HS spin chains of Sat is a further indication of the existence of a motif - 
based description of their spectrum, which would make possible a systematic study of 
the thermodynamics of these chains along the lines of Ref. [50] . 
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Figure 2: Probability density function histogram of the level density of an su(2) ferromagnetic HS chain 
of BCm type with parameter (3 = 2 for N = 12 spins (left) and N = 20 spins (right), compared to a 
normal distribution with the same mean and standard deviation as the spectrum (continuous red line). 
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